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ABSTRACT 

We present interferometric observations in the ^^CO (2-1) hne and at 1.3 mm 
dust continuum of the low-mass protostellar binary system in the cometary glob- 
ule CG30, using the Submillimeter Array. The dust continuum images resolve 
two compact sources (CG SON and CG SOS), with a linear separation of ~ 8700 AU 
and total gas masses of ~ 1.4 and ~ 0.6 Mq, respectively. With the CO images, we 
discover two high- velocity bipolar molecular outflows, driven by the two sources. 
The two outflows are nearly perpendicular to each other, showing a quadrupolar 
morphology. The northern bipolar outflow extends along the southeast (red- 
shifted, with a velocity up to ~2Skms~^) and northwest (blueshifted, velocity 
up to ~SOkms~^) directions, while the southern pair has an orientation from 
southwest (blueshifted, velocity up to ~lSkms~^) to northeast (redshifted, ve- 
locity up to ~41kms~^). The outflow mass of the northern pair, driven by the 
higher mass source CG SON, is ~ 9 times larger than that of the southern pair. 
The discovery of the quadrupolar molecular outflow in the CGSO protobinary 
system, as well as the presence of other quadrupolar outflows associated with 
binary systems, demonstrate that the disks in (wide) binary systems are not 
necessarily co-aligned after fragmentation. 

Subject headings: binaries: general — ISM: jets and outflows — ISM: individual 
(CGSO, BHR12) — stars: formation 



1. INTRODUCTION 

Binarity/multiplicity is frequent among protostellar systems (see e.g., Looney et al. 
2000). To probe these protostellar binary (protobinary) systems, high angular resolution 
observations of gas and optically thin dust emission at (sub-) millimeter (mm) wavelengths 
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are needed, because protostellar systems are deeply embedded in molecular cloud cores and 
are surrounded by infalling envelopes and accretion disks. However, these observations were 
long limited by the low angular resolution of single-dish telescopes and have only become 
possible with the recent availability of large (sub-) mm interferometers. During the past 
few years, interferometric mm observations have enabled us to discover several protobinary 
systems (Launhardt 2004), and there is an increasing number of interferometric studies of 
binarity in the protostellar stage (e.g., Chen et al. 2007). Nevertheless, due to the large effort 
and observing time required, the number of known and well-studied protobinary systems is 
still very small. 

CG30 (also known as BHR12, Bourke et al. 1995; or DC 253.3-1.6, Hartley et al. 
1986) is a bright-rimmed cometary globule located in the Gum Nebula region at a distance 
of ~400pc (Reipurth 1983). With high resolution submm (SCUBA; Henning et al. 2001) 
and mm (ATCA; Chen et al. 2008, hereafter Paper I) dust continuum observations, we have 
identified in this globule a wide protobinary system with a projected separation of ~21'.'7 
(~8700AU) and bolometric luminosities of 13.6 ±0.8 and 4.3 ± 0.5 for the northern and 
southern protostars, respectively. The two protostars appear to drive their own bipolar 
protostellar jets, which are roughly perpendicular to each other, as seen in the H2 1-0 5'(1) 
(Hodapp & Ladd 1995) and the Spitzer IRAC 4.5 /im images (Paper I). However, information 
on the molecular outflows in this protobinary system has been limited, although single-dish 
observations in the CO line (Nielsen et al. 1998) and CS line (Garay et al. 2002) had 
shown a large-scale (unresolved) bipolar outflow in the east-west direction. In this Letter, 
we report Submillimeter Arrajo (SMA; Ho et al. 2004) observations in the ^^CO (2-1) line 
towards CG 30 with angular resolution sufficient to allow us, for the first time, to map the 
distribution of outfiowing molecular gas from each member of the protobinary system. 



2. OBSERVATIONS AND DATA REDUCTION 

Observations at the SMA were carried out on 2008 March 7th in the compact configura- 
tion. Eight antennas were used in the array, providing baselines with projected lengths from 
14.1 to 136 m. The digital correlator was set up to cover the frequency ranges 219.5—221.3 
GHz and 229.5—231.3 GHz in the lower and upper sidebands, respectively. This setup in- 
cludes three isotopic CO lines of ^^CO (2-1) (230.538 GHz), ^^CO (2-1) (220.399 GHz), and 
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C^^O (2-1) (219.560 GHz), which were observed with an uniform channel spacing equivalent 
to 0.5kms~^. The 1.3mm dust continuum emission was also recorded with a total band- 
width of ~ 3.5 GHz (~ 1.7 GHz USB and ~ 1.8 GHz LSB). System temperatures ranged from 
200 to 400 K (depending on elevation), with a typical value of ~ 300 K. The primary beam is 
about 55" at 230 GHz. The quasar 3C84 was used for bandpass calibration, and the quasars 
0747-331 and 0727-115 for gain calibration. Titan was used for absolute flux calibration, 
which we judge accurate to ~20%, by comparison of the final quasar fluxes with the SMA 
calibration database. 

The data were calibrated using the IDL MIR package (Qi 2005) and imaged using 
the Miriad toolbox (Sault et al. 1995). With robust uv weighting 2, cleaned maps of the 
^^CO (2-1) emission and the continuum emission were produced with effective synthesized 
beam sizes of 5'.'5 x 2'.'4 and 5'.'0 x 2'.'0, respectively. The la rms noise levels are ~0.13 
(width ~0.5kms~^) and ~ 0.10 Jybeam"^ (width ~2.0kms~^) per velocity channel, and 
~4mJybeam~^ for the continuum emission. The final images were uncorrected for the 
primary beam attenuation. Further analysis was performed with the GILDA^ software 
package. 



3. RESULTS 

The SMA 1.3 mm dust continuum image of CG30 (Fig. 1) shows two compact sources, 
which were previously detected in our ATCA 3 mm dust continuum image (see Paper I). The 
angular separation between the two sources is measured to be 21'.'8 ± 0'.'5, consistent with the 
ATCA result. Following Paper I, we refer to the northern source as CG30N (also consistent 
with IRAS 08076— 3556) and to the southern source as CG30S. From Gaussian uv plane 
fitting, we derive flux densities of 256±51mJy for source N and of 99±20mJy for source 
S. The large-scale common envelope (see Fig. 1), detected in the single-dish 1.3mm maps 
with a total flux of ~ 1500 mJy (R. Launhardt et al., in preparation), is resolved out by the 
interferometer, namely that ~ 75% flux is missing in the SMA map. 

From the continuum fluxes measured at A 1.3 mm and 3 mm, we derive a spectral index 
a ~ 3.2 (where flux oc z/") in both sources, suggesting optically thin emission. The total gas 
mass in the circumstellar envelope was then calculated with the same method described in 
Launhardt & Henning (1997). In the calculations, we adopt an interstellar hydrogen-to-dust 
mass ratio of 110 (Draine & Lee 1984) and a factor of 1.36 accounting for helium and heavier 
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elements. For the two sources, we use a mass-averaged dust temperaturqjof Td = 20 K and an 
opacity of = 0.8cm^g~\ a typical value suggested by Ossenkopf & Henning (1994) 

for coagulated grains in protostellar cores. From the 1.3 mm dust continuum image, the 
total gas masses of sources N and S are estimated to be 1.4 ± 0.3 and 0.6 ± 0.1 M©, 
respectively, for the assumed distance of 400 pc. By contrast with the gas masses derived 
from the 3 mm dust continuum images, we estimate an opacity index of /? ~ 1.5 {ki, cx u^) for 
the two sources, a fairly typical index found in protostellar cores (see Ossenkopf & Henning 
1994). 

Figure 2 shows the velocity channel maps of the ^^CO (2-1) emission. For display pur- 
pose, we have smoothed the channel maps into a velocity resolution of 2.0kms~^. The 
CO emission is detected from 14jSR = — 26kms~^ to 48kms~^, with systemic velocity being 
~6.5kms~^. We note that the emission near the systemic velocity is scattered over the field 
of view, because the extended emission has been resolved out by the interferometer. Also 
note that the emission in the northeast and southeast corners of the panels between 8 and 
18kms~^ is aliasing due to insufficient uv coverage. In each panel, the two crosses indicate 
the positions of the two dust continuum sources. The CO emission from the two sources can 
be clearly distinguished due to their wide separation. For source CG 30N, blueshifted emis- 
sion (from Vlsr = ~ — 24kms~^ to ~3kms~^) is extending to northwest, while redshifted 
emission (from ~8kms~^ to ~29kms~^) is to southeast; For source CG30S, blueshifted 
emission (from ~ — 6kms~^ to ~5kms^^) is extending to southwest, while redshifted emis- 
sion (from ~30kms~^ is to ~48kms~^) to northeast. 

Figure 3 shows the integrated intensity maps of the CO (2-1) emission, plotted on the 
Spitzer IRAC band 2 (4.5 /im) image. Two separated bipolar outflows, with projected 
lengths of ~ 27,000 AU (northern pair) and ~ 20,000 AU (southern pair), are clearly seen in 
the images. The position angles (P. A., measured with respect to the red lobe) are ~ 128±5° 
for the northern pair and ~57±3° for the southern pair, i.e., the two bipolar outflows are 
nearly perpendicular to each other. In the northern bipolar outflow, the blueshifted lobe is 
leaf-shaped and spatially coincident with the Herbig-Haro (HH) object HH 120 (see Fig 3b), 
while the redshifted lobe is double-peaked and pointing to the H2 knot No. 6 (Hodapp & 
Ladd 1995). H2 knot No. 2, located to the northwest of the blueshifted lobe, is also likely to 
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one reference temperature for both sources. The actual dust masses must be scaled with the (not well-known) 
dust temperature. 
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be associated with the CG 30N outflow. Centered at CG 30S is a highly coUimated bipolar 
jet seen in the IRAC 4.5 ;um image. The bipolar CO outflow driven by CG30S is spatially 
coincident with this jet and also shows a coUimated morphology, in particular the redshifted 
lobe (length to width ratio is ~4). Its blueshifted lobe is spatially coincident with the near- 
infrared emission (knot No. 3), with the dust continuum source being located at the apex 
of both the CO and infrared emissions (Fig. 3c). This is consistent with our suggestion in 
Paper I that the infrared emission at CG 30S is due to scattered hght in a cavity evacuated 
by the jet/outflow. 



4. DISCUSSION 

4.1. Outflow Physical Parameters 

Thanks to the versatility of the SMA correlator, the ^^CO (2-1) line was observed simul- 
taneously with the ^^CO (2-1) hne towards CG 30, which allows us to estimate the outflow 
mass. For source CG 30N, the ^^CO (2-1) blueshifted emission is detected at velocities from 
~ —2 to ~3kms~^, while redshifted emission is from ~8 to ~ 12kms~^; For source CG30S, 
however, no ^^CO (2-1) outflow gas is found, although weak emission is detected around sys- 
temic velocity. The outflow masses are then calculated following Cabrit & Bertout (1990): 
for the southern bipolar outflow, as well as high-velocity component of the northern pair 
where ^^CO (2-1) is not detected, the ^^CO (2-1) emission is assumed to be optically thin, 
and the mass is calculated at each channel where outflow emission is detected, and then 
summed; for the low-velocity component of the northern pair, a correction for line opacity 
is performed at velocities where both the ^^CO/^^CO (2-1) emissions are detected. In the 
calculations, we assume LTE conditions, an excitation temperature of 20 K, an abundance 
ratio of [H2]/[^'^C0] equal to 8.9 x 10^ (Cabrit & Bertout 1992), and an isotopic ratio of 
[^^CO]/[^^CO] = 71 (Wilson & Rood 1994). No correction for inclination effects was applied 
in the calculations. 

The derived outflow masses, as well as other outflow properties (i.e., momentum P, 
energy E, mass- loss rate Mout) force F^, and mechanical luminosity L^), are listed in Ta- 
ble 2. A simple comparison between the two bipolar outflows shows that the northern pair, 
driven by the higher (circumstellar envelope) mass source CG 30N, is much stronger than the 
southern pair (~9 times in mass and ~3 times in energy). The outflow parameters derived 
in CG30 are relatively small in comparison with the general range of such parameters for 
sources of similar mass (see Wu et al. 2004). Nevertheless, we want to mention that the 
parameters in Table 2 only refer to the compact outflows detected in the SMA maps. The 
nature of the "extended outflow" components is unclear yet and requires confirmation by 
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combining the SMA data with single-dish data to recover extended missing flux. Therefore, 
these parameters represent lower limits. Furthermore, as suggested by the two larger-scale 
jets, mosaicked CO observations are probably needed to fully map the two outflows. 

4.2. Binarity v.s. Quadrupolar Outflows 

The two bipolar outflows found in CG 30 resemble a combined quadrupolar morphology. 
A similar outflow morphology was also found in several other (low-mass) sources, like e.g., 
IRAS 16293-2422 (Mizuno et al. 1990), HH288 (Gueth et al. 2001), and L723 (Lee et al. 
2002 and references therein). It is important to note that in all these sources protobinary 
systems have been discovered by high-resolution interferometric observations (IRAS 16293, 
Mundy et al. 1992; HH288, Gueth et al. 2001; L723, Launhardt 2004). This strongly 
supports the hypothesis that many quadrupolar outflow^ are actually two independent 
outflows driven by a binary system (Anglada et al. 1991), although several other scenarios 
have been proposed (see Arce et al. 2007). Furthermore, assuming that binary systems 
are formed through the fragmentation of collapsing protostellar cores (see Goodwin et al. 
2007), these quadrupolar outflows demonstrate that the disks in (wide) binary systems are 
not necessarily co-aligned after fragmentation. 

On the other hand, if binarity among protostars occurs as frequently as among pre- 
main sequence stars, one would expect to find many quadrupolar /multipolar protostellar 
outflows. Such kind of outflows are however rare in our observations. One possible reason is 
that protostellar disks (jets) are preferentially aligned during fragmentation, and they may 
then drive one common large-scale outflow. An example is the source HH 30 (see Guilloteau 
et al. 2008 and references therein). At the same time, the outflow strength of protostars 
is proportional to the circumstellar envelope mass (Bontemps et al. 1996). Since binary 
protostars appear to have very unequal circumstellar masses (Launhardt 2004), their relative 
outflow strengths are expected to be also very unequal. Therefore, another possible reason 
is that the weaker outflow may often not be detected with current single-dish telescopes that 
do not have sufficient angular resolution and sensitivity. The discovery of an equal-mass 
protobinary system in L 723 (Launhardt 2004) and IRAS 16293 (Mundy et al. 1992) and of 
a weak secondary outflow from the lower-mass component in the unequal-mass protobinary 
systems CG30 (this work), CB230 (Launhardt 2004), and BHR71 (Bourke 2001, Parise et 
al. 2006) supports this explanation. However, high-resolution outflow data so far exist only 
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for very few systems, and it is unclear yet why multipolar protostellar outflows are so rare 
in our observations. 

We thank the SMA staff for technical support during the observations, and Mark Gur- 
well for his maintenance of the Submillimeter Calibration List. Support for this work was 
provided by NASA through contracts 1224608 and 1279198 (T. L. B.) issued by JPL, Caltech, 
under NASA contract 1407. 
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Table 1. SMA 1.3 mm dust continuum results of CG30 



Source 




FWHM 


sizes" 






[mJy] 


maj . X min. 


P.A. 


[Mq] 


CG SON 


256±51 


3'.'4xiy8 


-86±5° 


1.4±0.3 


CG 308 


99±20 


5'.'5x2'.'3 


11±7° 


0.6±0.1 



"Flux densities and FWHM sizes of the continuum 

sources derived from Gaussian uv plane fitting. The er- 
ror bar of flux density is derived from yj c^^jj + o-g^ , where 
fcali is the uncertainty from calibration (~ 20% of flux den- 
sity) and fTfit is the uncertainty from Gaussian fitting. 

''Total gas mass; See text for dust temperature and opac- 
ity used. 



Tabic 2. Outflow paranw^tn-s 





Mass" 


Momentum" 


Energy" 


Mont'' 


Force*" 


Luminosity'' 


Gomponent 


[10-3 M0] 


[Mglcms-l] 


[10'*3 ergs] 


[lO-6M0yr-l] 


[lO-SM0kms-Vr-^] 


[Le] 


CG30N blue 


17.4 


0.19 


3.1 


7.3 


8.1 


0.11 


red 


37.9 


0.48 


8.7 


16 


19.8 


0.30 


CG 30S blue 


3.6 


0.047 


0.6 


2.1 


2.7 


0.03 


red 


2.3 


0.095 


4.0 


1.3 


5.6 


0.19 



''The lower limits derived from the SMA observations. 



The outflow mass-loss rate, force, and mechanical luminosity are estimated from the mass, momentum, and energy 
with dynamical ages (~2400yr for the northern pair and ~1700yr for the southern pair). These dynamical ages only 
refer to the compact outflows detected in the SMA maps. 
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8''09™34^0 33^0 32"0 3l"0 

R.A. (J2000) 



Fig. 1. — SMA 1.3mm dust continuum image of CG30. Contours start at ^3 a 
(l(T~4mJybeam~^) with steps of ~2(T to 17 cr, and then increase with steps of ~5(T. 
The grey dashed contour represents the half-maximum level of the 1.3 mm continuum emis- 
sion observed with SEST (Launhardt et al. in prep.). The synthesized SMA beam is shown 
as a grey oval in the bottom right corner. 
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Fig. 2. — Velocity channel maps of the ^^CO (2-1) emission (phase center R.A. =08:09:32. 998, 
DEC=-36:05:08.00, J2000). The channels were smoothed into 2.0 km s"^ per step from the 
original velocity resolution of 0.5kms~^. The center velocity of each channel is written in 
top left corner of each panel (in kms~^). The systemic velocity of the molecular cloud is 
~6.5kms~^. Contours levels correspond to 3, 6, 10, 15, and 20(7, then increase in steps of 
10(7, where the la level is ~0.1Jybeam~^. In each panel, the crosses mark the positions of 
the two dust continuum sources, the open dotted circle shows the primary beam (~ 55" at 
230 GHz), and the filled ellipse (lower right corner) indicates the synthesized beam of SMA. 




Fig. 3. — (a) The integrated intensity maps of the ^^CO (2-1) emission, plotted on the 
Spitzer IRAC2 (4.5;Ltm) image of CG30 (from Paper I). The SMA CO contours (color) 
are described below. The dotted circle shows the primary beam of SMA. Sources CG30N 
and CG 30S are labeled as "N" and "S" , respectively. The two dashed arrows show the 
directions of the infrared jets, while the arrow delineating the CG30N jet, redefined by the 
bipolar CO outflow, is ~ 20 degrees different to the one marking the same jet shown in 
Fig. 7 of Paper I. The knots in the jets are labeled with the same numbers as in Hodapp 
& Ladd (1995); (b) An enlarged view of the IRAC2 image, overlaid with the SMA 1.3 mm 
dust continuum contours (white) and the ^^CO (2-1) emission contours (reference position 
at R.A.=08:09:04.082, DEC=-36:00:53.53, J2000). The map shows the bipolar outflow 
driven by CG 30N. The blue (red) contours represent emission integrated over the velocity 
range — 24kms~^ <VLSR<3kms~^ (8kms~^ < Vlsr < 29 kms~^), which is blueshifted (red- 
shifted) respect to the cloud systemic velocity. SMA CO contours are from 10 Jykms"^ by 
step of lOJykms"^. (c) Same as Fig. 3b, but the map shows the bipolar outflow driven 
by CG30S. The blue (red) contours represent emission integrated over the velocity range 
— 6.5kms~^ < Vlsr < 5kms~^ (30kms~^ < Vlsr <48kms~^). SMA CO contours are from 
5 Jykms^^ by step of 5 Jykms"^. For display purpose, we use yellow contours in Fig. 3a to 
show the redshifted emission from source CG 30S. 



